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摘 要 : 气孔 时 度 (g,) 是 衡量 植物 与 外 界 环境 中 水 .CO; 等 物质 交换 速率 的 重要 参数 ,其 观测 与 模拟 
可 以 有 效 指示 物质 交换 情况 及 各 项 生理 参数 。 通 过 LI-6400 人 和 便携 式 光 合 测 定 仪 对 黄土 高 原 苹 果树 
叶片 各 项 生理 参数 进行 观测 ,分 析 g, 日 变化 特征 及 其 与 环境 因子 的 响应 关系 ,运用 Jarvis 模型 和 
Ball-Woodrow-Berry(BWB ) 模 型 对 g, 进 行 模拟 。 研 究 结果 显示 :(1) 黄土 高 原 苹果 树 &, 日 变化 在 气 
温 较 高 、 辐 射 较 强 的 8、9 月 观测 日 内 呈 双 峰 曲 线 。 上 午 (8:00 一 12:00) 随 着 太阳 辐射 逐渐 增强 ,气孔 
张 开 ,g, 在 11:00 一 13:00 出 现 第 一 次 峰值 ;中 午 (12:00 一 14:00) 由 于 气温 (ZT,) 升 高 ,为 避免 细胞 散失 
过 多 水 分 ,气孔 关闭 出 现 短暂 的 “光合 午休 ”现象 ;午后 (14:00 一 18:00), 随 着 7T,、 光 合 有 效 辐射 
(PAR) 下 降 ,g, 逐 渐 增加 并 在 15:00 一 17:00 出 现 第 二 次 峰值 (2) 通过 灰色 关联 度 分 析 ,发 现 g, 与 各 
环境 因子 的 关联 程度 依次 为 :PAR(0.731)>C0; 浓 度 (C,,0.712)> 饱 和 水 汽 压 差 (VPD,0.702)>7, 
(0.689)> 相 对 湿度 (有 ,0.673)。g. 与 各 环境 因子 的 响应 关系 表现 为 : 随 PAR、T,、C,h, 的 增 大 而 增 大 ， 
VPD 的 增 大 而 减 小 ,(3) 从 g, 的 模拟 结果 可 以 看 出 ,Jarvis 模 型 的 决定 系数 (0.678)、 修 正 效 率 系数 
(0.335) 和 修正 一 致 系数 (0.803) 均 优 于 BWB 模型 各 值 (0.329、-1.630、0.138), 而 平均 绝对 误差 
(0.103) 小 于 BWB 模型 (0.143) ,表明 Jarvis 模型 模拟 效果 较 好 。 通 过 对 黄土 高 原 苹 果树 叶片 &, 环 境 
响应 与 模拟 的 分 析 研 究 ,对 于 掌握 苹果 树叶 片 一 天 内 不 同时 刻 对 水 分 的 需求 变化 ,进一步 提高 该 


地 区 苹果 树种 植 的 水 资源 利用 效率 具有 重要 意义 。 
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气孔 是 调控 植物 与 外 界 环 境 中 水 和 CO0; 等 物质 
能 量 交 换 的 关键 环节 ,其 变化 深刻 影响 着 植物 蒸腾 
和 光合 作用 ,并 对 空气 相对 湿度 Ch., % ) 和 气温 
(ZT,,"C) 产 生 重要 的 调节 作用 "。 人 气孔 交换 速率 通常 
用 气孔 导 度 (g,,mol*H,0.m*…s ') 衡 量 , 植 物体 通过 
& 实 现 对 光合 作用 .蒸腾 作用 等 生理 特性 的 影响 抑 
制 ,并 随 7,、CO; 浓 度 (C,, hmol*mol') 等 条 件 的 变化 
而 发 生 改 变 , 在 植物 体 中 起 着 重要 的 平衡 调节 作 
用 。 通 过 分 析 g 与 各 环境 因子 的 响应 关系 ,根据 植 
物 的 生理 需求 来 调控 适宜 的 生长 环境 是 提高 植物 
水 分 利用 效率 的 重要 途径 ,对 于 作物 水 分 利用 机 制 
过 程 的 解释 以 及 环境 的 优化 管理 等 都 具有 重要 意 


收 稿 日 期 : 2020-10-08; 修订 日 期 : 2020-11-23 


义 ”。 在 & 的 观测 过 程 中 ,由 于 所 需 工 作 量 大 且 容 
易 受 到 外 部 天 气 状 况 的 影响 ,导致 在 获取 长 时 间 序 
列 的 数据 方面 存在 困难 。 因 此 ,运用 模型 进行 模拟 
成 为 研究 植物 g. 最 有 效 的 方法 。 目 前 广泛 用 于 估算 
名 对 环境 因子 响应 的 经 验 模型 主要 为 Jarvis 模 型 和 
Ball-Woodrow-Berry (BWB ) 模 型 ? 。 

苹果 是 一 种 重要 的 经 济 作物 ,开展 苹果 树 & 与 
各 环境 因子 的 啊 应 及 模拟 研究 ,对 于 提高 该 地 区 苹 
果树 种 植 的 水 资源 利用 效率 、 进 行 g, 尺 度 扩 展 及 其 
与 大 气 间 的 通 量 估算 等 都 具有 重要 意义 。 有 关 蔷 
果树 &, 的 变化 趋势 与 环境 响应 研究 , 张 敏 “ 对 不 同 
灌水 量 . 蔓 水 坑 深 条 件 下 苹果 树 & 日 变化 特征 以 及 
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& SICA TE (P,, umol. m°. s~) AMA Ee FA EAT 
FURS OT, OC UE DEE Te KHER FE 
FERR o 5 25 EER (T., mom”. s) AY H Aet 
征 及 二 者 的 响应 关系 ; 李 波 ”\ 李 营 ” 通 过 分 析 & 与 
土壤 含水 率 .光照 强度 OGG A CLT (PAR, W:m?) 
FA) Mh NYG FA, BE — AE EE TAS Te SER AB THE KE 
PRR RLF Fg. EA SAS EK 5 OCIA SOT eK 
SUE ARE FERRE g HEERE ETT T 
研究 ,结果 均 表 明 苹 果树 g, 日 变化 呈 先 升 高 后 降低 
的 趋势 ,个 别 观测 日 内 有 “光合 午休 ”现象 。 在 g. 的 模 
拟 研 究 方 面 , 郑 利 剑 ” 运 用 Jarvis 模 型 对 蓄 水 坑 灌 条 
EB ERY g 进行 了 数值 模拟 ,高 照 全 等 … 利 用 
Leuning(1995 ) 改 进 的 g 半 机 理 模型 模拟 了 不 同治 
水 条 件 下 侈 裁 苹 果树 7、g,、P, 的 动态 变化 。 关 于 苹 
果树 g 的 模拟 均 是 基于 单一 模型 ,缺乏 多 模型 模拟 
精度 对 比 研 究 。 

黄土 高 原 属于 典型 的 温带 季风 气候 区 ,区 内 降 
水 集中 、 四 季 分 明 , 是 苹果 树种 植 的 优势 产 区 。 然 
而 ,该 区 冬 春季 节气 候 干 旱 , 降 水 分 布 不 均 导 致 苹 
果树 生长 所 需 的 水 分 、 养 料 等 供应 不 足 ,使 其 生长 
受到 限制 。 本 文 以 黄土 高 原 苹果 树 为 例 , 研 究 苹 果 
树 & 与 各 环境 因子 的 啊 应 关系 ,选用 Jarvis 模 型 和 
BWB 模 型 对 黄土 高 原 苹 果树 &, 进 行 数值 模拟 , 旨 在 
擎 握 黄 土 高 原 苹 果树 一 天 内 不 同时 刻 对 水 分 的 需 
求 变化 ,进一步 提高 该 地 区 苹果 树种 植 的 水 资源 利 
用 效率 。 


1 研究 区 概况 


试验 样 地 位 于 山西 省 农 科 院 植 物 保护 研究 所 
东阳 试验 研究 基地 ,地 处 山西 省 中 部 地 形 区 ,属于 
温带 季风 气候 区 ,夏秋 季 降 水 集中 , 冬 春 季节 较 干 
FA ,平均 年 降水 量 为 450 mm, hk 60% ,年 均 气温 为 
9.5 ,年 均 风 速 为 2.4 ms'!。 试验 样 地 内 共有 苹果 
树 42 棵 ,平均 冠 幅 为 316 cmx318 cm ,平均 地 径 为 
10.9 cm。 


2 研究 方法 


2.1 试验 设计 与 方法 

试验 于 2019 年 7 一 9 月 苹果 树 主要 生长 季 内 进 
行 , 在 试验 样 地 内 选取 3 株 长 势 较 好 的 苹果 树 作 为 
样 树 ,采用 美国 LI-6400 便携 式 光 合作 用 测定 系统 


对 P,、T、C、 叶 表面 C0; 浓度 (C,, mol* mol') .PAR ,g, 
等 生理 参数 进行 观测 。 观 测 过 程 中 ,在 样 树 上 选取 
3 片 阳 生 叶 片 并 对 各 参数 值 重复 记录 3 次 ,观测 时 间 
为 8:00 一 18:00 ,时 间 间 隔 为 1ph, 不 同 观测 日 期 再 根 
据 天 气 .日 出 日 落 时 间 等 变化 情况 做 出 相应 调整 。 
此 外 ,为 获取 样 地 内 饱和 水 汽 压 差 (VPD ,kPa) 了 sh, 
和 风速 (wm's-) 等 气象 参数 ,还 设置 了 小 型 气象 站 
以 供 观测 。 数 据 收 集 整理 完成 后 ,采用 SPSS 24.0 
对 数据 进行 预 处 理 ,方便 后 期 利用 模型 进行 模拟 
分 析 。 
2.2 g; 模 型 介绍 

Jarvis 模型 综合 考虑 了 PAR、VPD、7, 和 C, 对 
植物 叶片 g, 的 影响 ,其 表达 式 如 下 : 

£,=8,(PAR)g,(VPD)g,(7,)g.(C,) (1) 

式 中 :g,(PAR)、gs(VPD)、gx(7T,)、g.(C) 分 别 表示 PAR 、 
VPD、T,、C, 对 g, 的 响应 函数 。 在 以 往 的 研究 中 , 国 
内 外 学 者 关于 单一 因子 对 g, 的 响应 模型 给 出 了 不 同 
的 表达 式 : 


gn(PAR)=a + a,PAR (2) 

_ aPAR 
gnaPAR)= PAR 3) 

_ a,t+a,PAR/a, 
&a(PAR)= Ty PARJa,_ (4) 
gm(D)=1 -026.(D) (5) 
1 

gm(D)= 1 +b,(D) (6) 

1-b,(D) 

D)= ! 7) 
gml ) 1+b,(D) ( 
grT) =c toT, +¢,T° (8) 

&(C,)=1-d4,C, (9) 


式 中 : an an a Dis bas cc cs d, BN IEE 
参数 ; D 为 VPD。 

BWB 模型 假设 g, 和 4 之 间 呈 线性 关系 ,在 很 
大 程度 上 描述 了 气孔 的 开 闭 机 理 , 计 算 公 式 如 下 
所 示 : 


Ah 
8,=m + go (10) 


式 中 :m 和 g, 为 待定 参数 ;4 为 光合 有 效 辐 射 ; C ,为 
叶 表 面 CO, 浓度。 
2.3 模型 模拟 精度 评定 指标 

通过 对 比 实测 叶片 气孔 导 度 与 各 模型 模拟 的 
气孔 导 度 结果 ,进而 可 以 确定 各 模型 的 模拟 精度 。 
本 文 引 用 Legates 等 “人 研究 中 所 选用 的 3 个 数据 指 
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WHE: p 为 分 辩 系 数 , 一 般 在 0~1 之 间 ,通常 取 0.5; 


一 致 系数 (d) 和 平均 绝对 误差 (MAE)。 各 系数 计算 
公式 为 : 


N 


2.0,- mM, 
E =1- #5 (11) 
$]o.-7; 
z 
20,- Mi 
d,=1-———=! (12) 
20, -0;|+|m,-0;) 
7 N 
2|0.-M 
MAE = ——— (13) 


式 中 :0 为 实测 叶片 气孔 导 度 (mol.HO.m2?.s 7) 5M, 
为 模拟 气孔 导 度 (mol.HO.m2.s ); 0, 为 实测 平均 
气孔 导 度 (mol.HO.m2.s);N 为 样本 总 数 。 如 果 
决定 系数 ( 尼 ) .El 和 di 的 值 较 大 而 MAE 的 值 较 小 ， 
且 和 斜率 接近 于 1, 则 该 模型 模拟 精度 高 。 
2.4 灰色 关联 度 

灰色 关联 度 分 析 主 要 是 根据 变量 间 变 化 发 展 
趋势 的 相似 或 一 致 性 程度 来 反映 变量 间 的 关联 程 
度 。 在 对 人 进行 灰色 关联 度 分 析 之 前 ,首先 要 建立 
一 个 灰色 系统 号 :选择 5$ 个 天 气 晴 朋 .数据 记录 较 好 
的 典型 日 作为 数据 源 , 将 g, 设 为 参考 序列 ,PAR、 
VPD、7T,、C, 和 及 这 5 个 环境 因子 设 为 比较 序列 ,进行 
关联 系数 (&;) 和 关联 度 (7;) 的 测定 。 
2.4.1 无 量 岗 化 处 理 ”由 于 影响 苹果 树叶 片 g, 的 各 
环境 因子 的 生理 意义 以 及 取 值 范围 等 都 有 较 大 的 
不 同 ,需要 对 数据 进行 无 量 纲 标 准 化 处 理 。 通 常 无 
量 纲 化 处 理 方法 包括 均值 化 法 , 初 值 化 法 和 xa 
变换 等 ,本文 主要 采用 变换 法 进行 数据 处 理 。 
2.4.2 关联 系数 和 关联 度 ”对 无 量 纲 化 后 的 数据 进 
行 灰 色 关联 度 分 析 , 通 过 计算 可 得 g, 与 各 环境 因子 
的 &; 和 x, 这 里 的 关联 程度 实质 上 是 指 曲线 趋势 发 
展 的 相似 程度 。 本 文选 定 g, 作 为 参考 序列 ,PAR、 
Ca Ta VPD .作为 比较 序列 ,计算 参考 序列 与 比较 
序列 在 各 个 时 刻 ( 即 曲线 上 各 点 ) 的 ,和 7,。 计 算 
方法 如 下 公式 : 


_ A(min) + pA(max) 
bn = A pA oa 
r= VEW (15) 


参考 序列 与 比较 序列 之 差 最 小 值 记 为 Amin) ,最 大 
值 记 为 Amax); 参 考 序列 与 比较 序列 每 一 点 的 绝对 
HEWA Ak), FER xi RR PAR .CT VPD 这 5 
个 比较 序列 ,表示 环境 因子 的 观测 点 (=1, 2, 
3, …) ;NN 为 各 环境 因子 的 观测 总 数 ;Ei( 有 为 比较 序 
列 中 各 时 间 点 的 值 。 


3 结果 与 分 析 


3.1 ZRH g AEE 

基于 2019 年 7 一 9 月 苹果 树 主 要 生长 季 内 实测 
& 数 据 , 每 月 选取 1~2 d 晴 朗 天 气 条 件 分 析 苹 果树 叶 
片 & 变 化 情况 (图 1)。 

由 图 1 可 以 看 出 ,苹果 树叶 片 & 日 变化 在 气温 
较 高 .辐射 较 强 的 8.9 月 观测 日 内 呈现 双 峰 曲线 。 
上 午 随 着 光照 的 不 断 增强 ,气孔 开 度 逐渐 增 大 并 在 
11:00 一 13:00 达 第 一 峰值 。 中 午 高 温 辐射 条 件 增 大 
了 叶 内 外 水 汽 压 梯度 (正午 气温 一 般 约 在 36.8 C, 
有 时 其 至 超过 40 % ,图 2) ,蒸腾 作用 也 逐渐 强盛 ,为 
减少 散失 过 多 水 分 ,气孔 闭合 ,g, 减 小 。 午 后 7 下 
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图 1 气孔 导 度 (g,) 日 变化 特征 


Fig. 1 Diurnal variations of stomatal conductance (g,) 
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图 2 气温 (Z) 日 变化 


Fig. 2 Diurnal variations of air temperature (T4) 
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降 ,g, 逐 渐 回 升 并 出 现 第 二 次 峰值 ,有 旦 下 午 g, 峰 值 始 
终 低 于 上 午 峰 值 。 

3.2 苹果 树 g, 的 环境 响应 

3.2.1 磋 色 关联 度 分 析 ”根据 无 量 纲 标 准 化 后 的 数 
据 分 别 计算 é; 和 x( 表 1 和 表 2)。 表 1 中 显示 了 各 个 
时 刻 g 与 PAR、C,、VPD、7T, 和 之 间 的 关联 程度 ， 
表 2 则 从 总 体 上 反映 g, 与 各 环境 因子 的 关联 度 。 


表 1 各 环境 因子 与 气孔 导 度 (g,) 的 关联 系数 
Tab.1 Correlation coefficient between each 


environmental factor and stomatal conductance (g,) 


PAR/W:m?* — C,/pmol+ mol" VPD/kPa TIC h, 1% 
0.732 0.423 0.533 0.665 0.484 
0.877 0.481 0.555 0.874 0.527 
0.815 0.524 0.818 0.677 0.693 
0.746 0.498 0.721 0.772 0.628 
0.838 0.407 0.751 0.714 0.516 
0.492 0.400 0.801 0.843 0.642 
0.998 0.522 0.513 0.510 0.590 
0.736 0.830 0.646 0.782 0.671 
0.911 0.969 0.958 0.775 0.897 
0.649 0.931 1.000 0.818 0.833 
0.864 0.804 0.893 0.882 0.723 
0.678 0.856 0.947 0.992 0.753 
0.834 0.521 0.476 0.419 0.497 
0.708 0.678 0.734 0.893 0.896 
0.612 0.611 0.887 0.915 0.698 
0.822 0.845 0.897 0.745 0.718 
0.668 0.850 0.702 0.537 0.729 
0.561 0.937 0.717 0.529 0.771 
0.854 0.646 0.817 0.810 0.645 
0.859 0.671 0.841 0.834 0.671 
0.755 0.716 0.486 0.510 0.558 
0.748 0.797 0.445 0.508 0.544 
0.856 0.802 0.334 0.413 0.468 
0.480 0.850 0.509 0.604 0.499 
0.356 0.702 0.646 0.790 0.599 
0.704 0.732 0.574 0.520 0.978 
0.787 0.714 0.648 0.541 0.774 
0.767 0.840 0.673 0.580 0.662 
0.864 0.966 0.932 0.770 0.592 
0.341 0.831 0.594 0.445 0.946 


iE: PAR 为 光合 有 效 辐 射 ;C, 为 C0, 浓 度 ;VPD 为 叶片 饱含 水 汽 压 
差 ;7, 为 气温 ;有 为 空气 相对 湿度 。 下 同 。 


由 表 2 可 知 ,g, 与 各 环境 因子 的 关联 程度 依次 
为 :PAR (0.731)>C, (0.712)>VPD (0.702)>T, (0.689)>h, 
(0.673). X LAA PAR, Ca, VPD 与 g, 的 关系 较为 密 
切 , 对 气孔 开 度 的 影响 较 大 ,而 有 与 g, 的 关联 度 
最 小 。 
3.2.2 g; 对 各 环境 因子 的 响应 ”太阳 光照 是 植物 进 
行 光合 .蒸腾 作用 的 重要 条 件 ,也 是 影响 植物 叶片 g&， 
的 重要 因素 。 由 图 3a 中 可 以 看 出 ,g, 随 PAR 的 增 大 
而 逐渐 增 大 。 上 午 (8:00 一 12:00)g, 随 PAR 的 增 大 而 
逐渐 增 大 ,并 在 11:00 一 13:00 间 出 现 第 一 次 峰值 ; 
午间 (12:00 一 14:00) PAR 持续 增强 , 当 PAR 超过 
1200 W- m° Ja (AF E] PAR 均值 达 1385 Wem”) ,为 减 
少 过 多 的 水 分 散失 气孔 闭合 ,g, 也 随 之 下 降 ; 午 后 
(14:00 一 18:00) 随 着 PAR .也 下降, 叶 水 势 得 到 一 定 
程度 的 回升 ,g, 也 呈 短 暂 的 上 升 趋势 并 出 现 第 二 次 
峰值 ,之 后 g, 趋 于 下 降 。 

C0; 是 植物 进行 光合 作用 的 重要 原料 ,通常 认 
为 低 C, 可 以 促进 气孔 开张 ,高 C, 会 使 气孔 迅速 关 
闭路 。 图 3b 给 出 了 g, 对 CC, 的 响应 关系 ,可 以 看 出 ， 
当 C, 较 低 时 ,g, 随 C, 的 增 大 逐渐 增 大 ; 当 C, 增 大 到 
375 umol.mol JA , g, BE C. 下降 的 趋势 较为 明显 。 维 
持 较 高 & 的 C, 范 围 约 在 350~375 jumol + mol! Z [i , 
过 高 或 过 低 的 C, 都 会 对 植物 g. 产 生 一 定 抑制 作用 。 

图 3c 给 出 了 g, 对 VPD 的 响应 关系 ,可 以 看 到 ,g， 
随 VPD 的 增 大 而 逐渐 降低 。 从 全 天 来 看 ,VPD 从 早 
间 开 始 持 续 上 升 ,g, 下 降幅 度 较 小 ;午间 VPD 迅速 增 
大 (午间 VPD 均 值 为 3.784 kPa, 峰 值 可 达 7.107 kPa), 
气孔 闭合 以 避免 过 多 的 水 分 散失 ,g, 值 也 迅速 下 降 。 

一 般 认 为 ,7, 对 植物 叶片 g. 的 影响 是 随 着 7 的 
升 高 ,g, 有 增 大 的 趋势 , 当 7, 超 过 某 一 界限 时 ,g, 迅 
速 降 低 。 由 图 3d 可 以 看 出 ,g, 随 7 的 升 高 逐渐 增 
大 ,在 30 C~35 左右 达到 最 大 值 ,此 后 随 7 的 升 高 
名 逐渐 降低 。 

图 3e 给 出 了 有 对 玉 的 响应 关系 ,可 以 看 出 ,g&, 随 
.的 增 大 而 增 大 。 从 全 天 来 看 ,清晨 及 较 大 (均值 约 
为 43.048% ,峰值 可 达 52.191%),g, 随 着 的 增 大 而 
上 升 ;午后 (14:00 一 18:00) 随 着 7T, 和 PAR 的 不 断 升 
高 ,及 迅速 下 降 ,g, 受 高 温 影 响 也 旺 下 降 趋势 。 


表 2 各 环境 因子 与 气孔 导 度 (8,) 的 关联 度 


Tab.2 Correlation between environmental factors and stomatal conductance (g;) 


关联 度 PAR/AW.mz= C,/mol+ mol™ 


VPD/kPa TIC hJ% 


Ti 0.731 0.712 


0.702 0.689 0.673 
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注 :g 为 苹果 树叶 片 气 孔 导 度 ;PAR DIAA CA COHEN | VPD Jy tL AA TERE T. h A CHAE 


图 3 苹果 树叶 片 气孔 导 度 (g,) 与 各 环境 因子 的 关系 


Fig.3 Response relationship between stomatal conductance (g,) and environmental factors in apple tree leaves 


3.3 黄土 高 原 苹 果树 g, 模 拟 
3.3.1 Jarvis 模 型 BWB 模型 表达 式 的 确定 由 
Jarvis 模 型 中 各 单一 因子 对 g, 的 响应 模型 可 以 看 出 ， 
由 不 同 因子 计算 公式 组 合 确定 的 Jarvis 模型 表达 式 
有 9 种 形式 。 本 文 基于 实际 观测 的 苹果 树叶 片 g, 数 
据 , 对 比分 析 Jarvis 模 型 不 同 表达 式 的 模拟 结果 , 比 
较 尼 值 确定 最 优 模型 表达 式 。 此 外 ,通过 SPSS 24.0 
对 观测 数据 进行 拟 合 确定 待定 系数 值 ,最 终 确定 
Jarvis 模型 和 BWB 模型 的 表达 式 ,具体 结果 见 表 3 
和 表 4。 

基于 确定 的 Jarvis 模型 和 BWB 模型 表达 式 , 结 
合 各 参数 值 模拟 计算 苹果 树叶 片 g,, 再 根据 评定 模 
型 精度 指标 ,综合 分 析 模 型 的 模拟 精度 。 
3.3.2 gala Rate 苹果 树叶 片 主要 生长 季 内 
实测 g 与 Jarvis 模 型 BWB 模 型 模拟 g, 对 比 见 图 4 和 


表 3 2019 年 用 于 模拟 苹果 树叶 片 气孔 导 度 
的 Jarvis 模型 结构 确定 
Tab. 3 Determination of the Jarvis model expressions 
for estimating stomatal conductance of 


apple tree in 2019, respectively 


表达 式 R 系数 

& (PAR) gp(VPD)g7(Z,)g.(C,) 7 S 
gun(PAR)g (VPD)g7(T)g.(C,) 0.677 0.997 
*g, (PAR) gp VPD)g(T,)8.(C.,) 0.678 0.756 
ga(PAR)g» (VPD)g7(T,)g.(C,) - - 
gn(PAR)g (VPD)g(T )g.(C,) 0.628 -0.222 
& (PAR)g,(VPD)g (7) (C,) 0.617 -0.861 
ga(PAR)g» (VPD)g7(T)g.(C,) - - 

& (PAR) 8 (VPD); T,)g C.) 0.673 0.517 
& (PAR) gp(VPD)g,(7,)g.(C,) 0.677 -0.262 


注 :* 表 示 该 模型 为 最 优 叶片 气孔 导 度 模型 ;- 表 示 无 法 进行 拟 合 
计算 。 
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#4 Jarvis 模型 .BWB 模型 表达 式 
Tab. 4 The expressions of the Jarvis and BWB models 


模型 公式 表达 式 R 
Jarvis HE — g,(PAR)g(VPD)gi(T)g(C.) g,=(0.001 +6.885E 7PAR) x 1—SS3OVED «(128,936 + 11.4447, 一 0.15072)x(L+0.003CJ) 0.678 
AXh AXh 
BWB 模 型 ” g,=m C “+g, g,=0.001 x C “+ 0.036 0.329 


表 5。 从 网 4 中 可 以 看 到 ,通过 Jarvis 模型 计算 得 到 3 PAR AT, AY 82 A], FL K IF, g RR EFIE 
的 全 果树 g, 模 拟 值 在 散 点 图 中 的 分 布 较 BWB 模 型 ” 11:00 一 13:00 达 第 一 峰值 ;午间 由 于 升 高 ,强烈 蔽 
更 为 集中 ,线性 拟 合 方程 斜率 (0.756) 也 高 于 BWB 腾 作 用 下 气孔 关闭 ,g, 减 小 ;午后 T 减 弱 ,g. 回 升 并 在 
方程 线性 拟 合 斜率 (0.321)。 从 表 5 中 可 以 看 到 ， ”15:00 一 17:00 出 现 第 二 次 峰值 ,有 旦 下 午 g, 峰 值 始 终 
Jarvis 模 型 的 (0.335) 、d1(0.803) 和 尺 (0.678) 均 高 ” 低 于 上 午 峰值 。 
F BWB 模型 各 对 应 值 (-1.630,0.138 和 0.329), 而 在 与 环境 因子 的 响应 关系 中 , 随 着 PAR .C, A. 
Jarvis 模 型 的 MAE(0.103) 低 于 BWB 模 型 (0.143)。 的 增加 ,叶片 & 均 表现 为 上 升 趋势 , 当 超过 各 因子 
因而 ,在 本 研究 中 ,Jarvis 模 型 对 苹果 树叶 片 & 的 模 ”一 定 阐 值 后 ,g 早 下 降 趋势 。 这 一 结论 与 柴 梦 溢 中 、 
拟 精度 更 高 。 高 娟 "对 苹果 树叶 片 g 与 PAR、7,, 以 及 司 建华 等 
ee 对 胡杨 g. 与 C,、 王 海 珍 等 ”对 灰 胡 杨 g 与 h 的 响应 
4 讨论 结果 一 致 ,但 与 阮 成 江 等 2 对 黄土 丘陵 沙 球 g 与 C,、 
本 文采 用 LI_6400 光合 作用 测定 系统 对 便 果 树 。 唐 风 德 等 ”对 长 白山 阔 叶 红 松 林 叶 片 g, 与 Ch 的 
& 进行 观测 ,分 析 睛 朗 天 气 条 件 下 & 的 日 变化 特征 ， 响应 结 末 存在 差异 ,这 可 能 与 研究 区 域 地理 环 境 和 
结果 显示 苹果 树 & 日 变化 在 气温 较 高 . 辆 射 较 强 ”植物 种 类 有 关 。 对 于 & 与 YPD 的 响应 , 结 末 显示 g 
的 8.9 月 观测 日 内 呈 双 峰 曲线 ,这 一 结论 与 高 照 全 ” 随 VPD 的 增 大 而 减 小 ,这 一 结论 与 王 海 珍 等 ”对 灰 
等 对 不 同 水 分 条 件 下 苹果 树 . 李 明 起 等 "对 成 果 ”胡杨 、 周 丽 娜 等 中 对 林 下 参 叶片 g, 与 VPD 的 响应 结 
末期 苹果 树 的 g, 日 变化 结果 一 致 ,其 原因 在 于 :上 午 。 果 一 致 。g, 对 环境 因子 不 同 的 啊 应 关系 原因 在 于 : 


0.35 [ (a) Jarvis 模 型 0.25 F (b) BWB 模 型 


本 e 

a ta e 
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x 0.20 e ee at's e8 e z 0.15 ee set e e.° 
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图 4 2019 年 苹果 树 主 要 生长 季 内 实测 气孔 导 度 与 模拟 气孔 导 度 对 比 
Fig.4 Comparison between the measured and simulated stomatal conductance during 
the main growing season of apple tree in 2019, respectively 
HS 2019 FE RHE Et KS Jarvis, BWB 模型 模拟 精度 对 比 
Tab. 5 Comparison of accuracies between the Jarvis and BWB models in 2019, respectively 
模型 修正 效率 系数 (E)/mm 修正 一 致 系数 (di)/mm 平均 绝对 误差 (MAE)/mm R 

Jarvis 模 型 0.335 0.803 0.103 0.678 


BWB 模 型 -1.630 0.138 0.143 0.329 
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一 天 之 中 , 随 着 PAR 7.C. 的 逐渐 增强 ,气孔 开 度 增 
大 导致 叶片 P, 加 快 ,g, 呈 上 升 趋势 ; 当 各 环境 因子 超 
过 一 定 阅 值 后 ,叶片 在 外 界 高 温 、 较 高 VPD 的 刺激 
下 ,为 减少 细胞 散失 过 多 的 水 分 ,气孔 闭合 ,g 呈 下 
降 趋 势 。 而 较 高 的 及 多 对 应 于 清晨 .晚间 或 阴云 . 光 
照 较 弱 等 天 气 状况 , 相 较 于 晴朗 天 气 条 件 下 的 g, 数 
值 较 低 。 

本 文 运 用 Jarvis 模 型 和 BWB 模 型 对 苹果 树 g, 进 
行 模拟 ,结果 表明 Jarvis 模型 模拟 效果 较 好 ,这 一 结 
Ves ABA GI iz AA Jarvis 模型 模拟 苹果 树 & 的 结论 
一 致 。 此 外 ,高 冠 龙 等 对 极端 干旱 区 胡杨 林 以 及 
马 蓉 等 汪 对 新 疆 博 斯 腾 湖 北岸 芦苇 叶片 人 的 模拟 结 
果 也 同样 显示 ,Jarvis 模型 的 拟 合 效果 更 好 。 未 来 应 
加 强 干 旱 半 干旱 气候 条 件 下 多 种 作物 的 g, 模 拟 以 及 
于 旱 区 自然 植被 蒸 散 量 估算 .模拟 和 生态 需 水 量 
计算 ” ,进一步 探讨 Jarvis 模型 的 适用 性 。 


5 结论 


x~ 


本 文 主要 对 黄土 高 原 苹 果树 g, 进 行 系统 论述 ， 
首先 分 析 了 叶片 g, 日 变化 特征 ,解释 了 g 与 各 环境 
因子 的 响应 过 程 ,运用 Jarvis 模 型 和 BWB 模型 进行 
模拟 分 析 , 得 到 的 研究 结论 如 下 : 

(1) 黄土 高 原 苹果 树 g, 日 变化 在 气温 较 高 . 辐 
射 较 强 的 8.9 月 观测 日 内 呈 双 峰 曲 线 , 且 上 午 峰 值 

台 终 高 于 下 午 峰 值 。 其 原因 在 于 :上 午 受 PAR 和 也 
影响 ,气孔 张 开 ,g& 也 迅速 上 升 并 在 11:00 一 13:00 达 
第 一 峰值 ;午间 由 于 7 升 高 ,强烈 的 蒸腾 作用 下 气 
和 孔 关闭 ,g, 减 小 ;午后 7 减弱 ,g, 回 升 并 于 15:00 一 17:00 
出 现 第 二 次 峰值 , 且 下 午 & 峰值 始终 低 于 上 午 峰 值 。 

(2) &, 与 各 环境 因子 的 关联 程度 依次 为 :PAR > 
C, > VPD > T, > h, Hg. BAPAR 、T,、C,h, 的 增 大 而 增 
大 , 随 VPD 的 增 大 而 减 小 。g, 对 环境 因子 不 同 的 响 
应 关系 原因 在 于 :一 天 之 中 , 随 着 PAR、7,、C 的 逐渐 
增强 ,气孔 开 度 增 大 导致 叶片 P, 加 快 ,g, 量 上 升 趋 
势 ; 当 各 环境 因子 超过 一 定 姜 值 后 ,叶片 在 外 界 高 
温 、 较 高 VPD 的 刺激 下 ,为 减少 细胞 散失 过 多 的 水 
分 ,气孔 闭合 ,g, 呈 下 降 趋势 。 而 较 高 的 及 多 对 应 于 
清晨 .晚间 或 阴云 .光照 较 弱 等 天 气 状 况 , 相 较 于 晴 
朗 天 和 气 条 件 下 的 g, 数 值 较 低 。 

(3) g, 的 模拟 结果 表明 Jarvis 模型 模拟 效果 较 
好 。 在 未 来 ,还 应 加 强 干 旱 半 和 干旱 气候 条 件 下 多 种 
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作物 的 & 模 拟 研 究 以 及 干旱 区 自然 植被 蒸 散 量 佑 
算 模拟 和 生态 需 水 量 计 算 ,进一步 探讨 Jarvis 模 型 
的 适用 性 。 
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Environmental response and modeling of stomatal conductance 


of apple trees on the Loess Plateau 
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Abstract: In this study, stomatal conductance (g,) was used to measure the exchange rate of water, CO», and other 
substances between plants and the external environment. The observation and simulation of g, effectively 
indicated the exchange of substances and various physiological parameters. The Loess Plateau is typical of a 
temperate monsoon climatic zone in China and supports an agricultural area dominated by the cultivation of 
apple trees. In this study, an LI- 6400 portable photosynthesis system was used to observe the physiological 
parameters of apple trees in situ on the Loess Plateau, the diurnal variation characteristics of their g., and the 
relationships of these variables with environmental factors. The g. was simulated using the Jarvis model and the 
Ball- Woodrow- Berry (BWB) model. The results show that (1) the diurnal variation of g, in apple trees on the 
Loess Plateau showed a bimodal curve in August and September when temperatures were high and radiation was 
strong. The solar radiation increased gradually in the mornings (8:00—12:00), the stomata opened, and the first 
peak of g, appeared between 11:00 and 13:00. Around the noon hour (12:00—14:00), stimulated by the increase 
in temperature (T}), the stomata of plants closed for a short period during the “midday depression of photosynthesis” 
to reduce water loss from plant cells. In the afternoon (14:00—18:00), as T, and photosynthetically-active 
radiation (PAR) decreased, g, gradually increased, and a second peak appeared between 15:00 and 17:00. (2) 
Using the gray correlation degree, the correlation between g, and various environmental factors was as follows (in 
descending order): (PAR, 0.731)>CO, concentration (C, 0.712)>vapor pressure deficit (VPD, 0.702)>7. (0.689)> 
relative humidity (A,, 0.673). The response relationship between g, and the various environmental factors 
produced the following results: (1) it increased with increases in PAR, Ta, Ca, and h, and decreased with increases 
in VPD, and (2) The simulation of g, showed that the value of the determination coefficient (0.678), the modified 
coefficient of efficiency (0.335), and the modified index of agreement (0.803) were higher in the Jarvis model 
than in the BWB model (0.329, —1.630, 0.138), and the mean absolute error (0.103) was smaller than the error in 
the BWB model (0.143). Comparison of the simulation accuracy of multiple models showed that the Jarvis model 
had superior simulation accuracy. The results of the analysis of g, in response to environmental factors and its 
simulation in apple tree leaves on the Loess Plateau is important to understand how the demand for water by the 
leaves changes throughout the day. This knowledge can be used to improve the efficiency of water utilization and 
thus optimize the harvest. In the future, simulation studies of g, for a variety of crops grown in arid and semi-arid 
climatic conditions may benefit from the application of the Jarvis model in resource management. 


Key words: stomatal conductance; grey relational degree analysis; apple tree; Loess Plateau 


